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Abstract 



^ \ We report numerical detection of new type of localized structures in the frame 

^ ■ of Maj da-McLaughlin- Tabak (MMT) model adjusted for description of es- 



sentially nonlinear gravity waves on the surface of ideal deep water. These 
structures - quasibreathers, or oscillating quasisolitons - can be treated as 



■ groups of freak waves closely resembling experimentally observed " Three Sis- 

ters" wave packs on the ocean surface. The MMT model has quasisolitonic 



solutions. Unlike NLSE solitons, MMT quasisolitons are permanently back- 
ward radiating energy, but nevertheless do exist during thousands of car- 
rier wave periods. Quasisolitons of small amplitude are regular and stable, 
but large-amplitude ones demonstrate oscillations of amplitude and spectral 
^ \ shape. This effect can be explained by periodic formation of weak collapses, 

carrying out negligibly small amount of energy. We call oscillating quasisoli- 
tons " quasibreathers" . 
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1. Introduction 

Development of analytic theory of freak (or rogue) waves is one of the most 
interesting problems of hydrodynamics. In spite of recent progress in this area 
[ij] many important questions are not answered yet. Apparently, freak waves 
are the structures well localized in space, see FigJT] But behavior of freak 
waves in time in co-moving coordinate frames is not still explored. From the 
experimental viewpoint this is a hard question. It cannot be answered by a 
resting observer, for whom the freak wave is just a single event localized in 
time, see FigJH From the other hand, satellites move too fast to record the 
full " live story" of a freak wave. 

The standard model for description of freak waves in deep water is Non- 
linear Schrodinger Equation {NLSE). This equation has a plethora of exact 
solutions which often are associated with the freak waves on deep water. 
Some of these solutions are presented in 2[, more recent developments can 
be found in articles |3|-[6|. These solutions, however, presume existence of 
background monochromatic wave (condensate) and are connected to the sub- 
ject of our paper only indirectly. For this reason, we do not pursue a purpose 
to present here the detailed description of all solitonic solutions on the con- 
densate background, as well as completed and controversial history of their 
discovery. In this article we study the solitons on almost zero background. 

However, we should mention the remarkable NLSE solution found by 



Peregrine 



7(1 . This solution in the co-moving coordinate frame is an instan- 



ton, describing the single event - appearance and disappearance of the freak 
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waves group. Today we can speak about two alternative versions of the freak 
wave theory. The "instantonic" version assumes that the freak wave is a sin- 
gle event, localized in time. The "solitonic" version proposes that the freak 
wave are described by persistent solitons, probably oscillating in time. So 
far experimental data are too scarce to make a conclusion in favor of one of 
these theories. 

One should remember that NLSE is derived in the assumption that the 
wave train size, containing the freak waves, is much larger than characteristic 
wave length. Most of collected experimental data, however, show that in the 
real ocean this condition is not satisfied (see FigJH [2]) and that the NLSE 
is hardly applicable. 

A level of nonlinearity of quasi-monochromatic wave group is measured 



by the characteristic steepness /i ~ /ca 
amplitude). Our numerical experiments 



k is the wavenumber and a is the 
8| show that NLSE is applicable if 



/i ^ 0.07. According to our calculations, NLSE is not applicable if /i ~ 0.1. 



Recent numerical experiments 



18| show that this limit might be extended 



to ^ 0.15. However, for freak waves in the real sea fi ~ 0.3 4- 0.5 (see 
Appendix II). The NLSE is absolutely not applicable for description of that 
steep freak waves. 

It is also known from observations [ll] that a typical configuration of a 
freak wave group consists of three sequent waves - "Three Sisters". This 
group is too short to be described by NLSE. 

What are the alternatives to NLSE model? The most consistent ap- 
proach is the use of the exact Euler equations for description of the potential 
flow of the ideal fluid with free surface. Some advances in this direction are 
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already achieved js]- 10|. However, the study of more simple and less accu- 
rate models also could be very useful. In this article we present our result on 



numerical solution of well-known MMT [ll|] equation with the special choice 
of parameter a = 1/2, /3 = 3, A = making this model well adjusted for 
description of surface gravity waves. 

Our results mostly support the "solitonic" theory of freak waves. We 
started with initial data, corresponding to NLSE solitons and discovered 
formation of persistent quasisolitons existing for more than two thousands of 
wave periods. These quasisolitons slowly radiate energy in backward direc- 



tion. As was shown recently 



12j in the "model case" (a = 1/2, /3 = 0, A = 1) 



this effect plays the key role in formation of the wave turbulent spectrum, 
but in our case its influence is negligibly small. 

However, we discovered completely new effect. While quasisolitons of 
small steepness ^ 0.1) behave similar to NLSE solitons on zero back- 
ground, the quasisolitons of higher steepness demonstrate almost periodic 
oscillations of amplitude and spectral shape, periodically forming power-like 
tails in spectra. This effect can be explained by modulational instability in- 
side the quasisoliton. Development of this instability leads to formation of 
"weak" one-dimensional collapses, which deform the spectrum, but absorb 
negligibly small amount of energy. Thereafter we call oscillating quasisolitons 
" quasibreathers" . 
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2. Basic model 



The Maj da-McLaughlin- Tabak (MMT) equation (see [llj, (l^ and (isj) 
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where ^^{x, t) is the complex function and the fractional derivative is defined 

by 
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(2) 



has been attracting lately fare attention of nonlinear wave scientists. The rea- 
son is MMT equation incorporates several already known important cases, 
and also can be used as a "test-bed" for verification of the concepts like weak- 
turbulent waves spectra, localized structures and their co-existence 1^, 15 |. 
For a = and (3 = 0,2 Eq. ([!]) is completely integrable. If a = 2 and (3 = 0, 
it is the classical NLSE for focusing (A = —1) and defocusing (A = -|-1) 
cases: 

,dip 



dt dx"^ 
If a = 2 and = 2, transformation 
tive NLSE [igI: 

.d(f) 
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iip turns Eq.([T]) into the deriva- 



dt 

Through Fourier transform 

ipk 
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ip{x)e dx 



Eq.dl]) can be rewritten in the form 
where 

Suppose that in Eq. (jl]) Tkkik2k3 is a generic function satisfying the symmetry 
conditions 

Tkk\,k2k2 Tkik,k2k3 Tkki,k3k2 Tk2k3,kki (6) 

For matrix coefficient ([5]) conditions (j6]) are satisfied and Eq.(fT]) is a Hamil- 
tonian system 



.dipk 


SH 


' dt = 


^Vk 


H = 


J, 



I 

'^k?dk + - 

Obviously, the Hamiltonian if is a constant of motion. Other motion con- 
stants are wave action 

N = j \i;k?dk 

and wave momentum 

Another model of type (jlj), describing surface waves on deep water, is so- 
called "Zakharov equation'' 



13| . This equation is not heuristic like MMT, 



it was systematically derived from Euler equations and therefore is supposed 
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to be more accurate in corresponding context. In this equation T^k,eki,ek2,tk3 = 
e^Tfefcifeafcs is cumbersome homogeneous function of the third order. 

One should note that if Eq. ([1]) is apphed for description of gravity waves, 
the surface shape can be reconstructed by the formula (see Appendix I) 

ri{x,t) = ^ [ e''^\k\'/\ilJk + rk)dk (7) 



V2 

3. Solitons and quasisolitons 

Let us look for a solution of Eq. (jl]) in a form 

Mt) = e'^^'-'^'Uk (8) 

where Q and V are the constants. The function (pk should satisfy the non- 
linear integral equation 

, .jTi234(l)t(l)2hS{k + ki-k2-h)dkidk2dk3 

'^^ = ^ -n + kv-\k\- 

This equation has solutions if Q and V can be chosen such that the denomi- 
nator in Eq. (|9]) cannot be zero for real k. This might happen only if a > 1. 
Let's suppose now a < 1. One can see that in this case the denominator 
in Eq. always has zero, which is clear from FigJSl Let Q < 0, V > 0. 
Thus, any solution of type ([8]) has singularity at negative k. It means that 
strict soliton solution of type ([8]) does not exist. However, one can construct 
approximate solutions, such that 0^ in dHD is slow function of time. These 
app^xtaate .ov.g =oMio,., .ad^ajing ene.gy ,„ the backward direction 
are called quasisolitons after paper [Ti 

As it was recently shown in jl^^lSj. quasisolitons play the central role 



in wave turbulence in the frame of MMT model if a = |, /3 = and A = +1. 
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It was shown that in this case the backward radiation plays the central role 
in dynamics of quasisolitons. But we study only the case a = |, /3 = 3 and 
A = +1. In this case , which intentionally models the gravity waves on deep 
water, the backward radiation is not that strong, due to essential nonlinearity 
suppression in the area of small wave numbers. Nevertheless, we definitely 
detect this phenomenon in our numerical experiments. 

Consider the structure of the denominator in Eq.Q. One can expect 
existence of the quasisoliton in the case when the straight line u = kV — Q 
is tangential to the curve u = k°'. The conditions of equal derivatives and 
existence of the common point of these two curves 8bt) k — k^yy, ctrci 

V = ak'^-' (10) 
Q = (a-l)C (11) 

We are now returning back to non-stationary Eq.(jl]) and make the change 
of variables k = k^ + k, k << k. Dispersion relation expansion into Taylor 
series 

(k^ + kT = kl + akl~'K + \a{a - l)k^-\' 
and change of variables 

^fc(t) = e-'(^™+"*^™"'")V«(^) (12) 



gives 



dt 2 



^a(a- 1)C~'/^V«+ (13) 



+ k^ j (t>*^^(t>^^(t)^^6{K + Ki- K2- K3)dKidK2dK3 = (14) 
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Another change of variables 



gives 



^%^ = ^«(a-l)C-'(g2 + «:^)x.+ (15) 

+ ik^ j XkiXk2Xk35{k + Ki- K2- K3)dKidK2dK3 

Applying inverse Fourier transform t) = J XK{t)e'^'^^dK, to the last equa- 
tion, we get NLSE in real space: 

+ - - - k'JxW = (16) 

Eq.( !T6|) has partial stationary solution 



which produces approximate quasisoliton solution of the Eq.([T]) with A = 1: 

^[x,t) = x{x - vt)e'^^+^^*e'''"^^''-'"^ (18) 

n = -(i-«)c 

A = -l«(l-a)C-V 



a-l 



V = ak,^ 

The characteristic wave-number k^ = —ck^ of backward radiation asso- 
ciated with the quasisoliton (see Figj3]) can be found from the equation 

koV-n=\ko\" (19) 
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together with Eq.([in])-(III]). For a = 1/2 

c = 3-V8~ 0.172 (20) 

Therefore, due to the smallness of the ratio -?";(/=o.fco,fcm.fcm) ^ ^3 ^ 5 . iq-3 
the backward radiation process in framework of the MMT model for /3 = 3 
is suppressed with respect to the case /3 = 0. 

To obtain the surface shape we replace in Eq.([7]) fc^^^ with kU^ and get 

^ = T TT" IT cos(wt - kmx) 

km cosnq[x — vt) 

Thus q is the standard steepness. 

4. Self-similar collapses 

Eq.(Il]) has self-similar solution: 

^(x,t) = (to-t)^/^F(^^^-^) (21) 
For the shape of the surface it gives 



At t —7- to time must vanish from Eq. ([T]), which means that 

7] — 7- a^x for a; > 
1] — )■ a~x for X < 

where a"*' > and < are the constants. In other words, solution Eq. 
( I2T]) describes formation of a wedge, in general (if a"*" 7^ a~), tilted with 
respect to vertical line. In fc-space we get 

^{k,t) = {to-tf'F{kito-ty) (23) 
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According to (1251) F{^) — )• ^ at ^ — )■ 0. Hence, asymptotically 

^/'(A;,t) ~ A;^^/^ (24) 
|^(A;,t)|2 ~ A;-9/2 ^25) 

Formation of collapses like Eq. fl2T]) -f l22]) means growth of power-like tails in 
A;-space. The spectrum Eq. fl2^ - fl2^ appears only at the moment of collapse 
t — i- 0. The singularity has the form of appearing and vanishing wedge, 
absorbing some amount of energy. However, time-averaged spectrum can 
have a slope different from IV^^P ~ A;"^/^. If the collapse events are rare, the 
slope must be higher that A;~^/^. 

Let us suppose that the collapse is "weak" and that only a very small 
part if energy is dissipated in an individual event. It means that the collapse 
is "almost" invertible process, symmetric in time with respect to the sign 
change to —t. In other words, the collapsing solution is 

ij{k,t) = \to-t\'/'F{k{to-t)') 

Now we can perform the Fourier transform in time and get 

muj) = |to - tf'F{k{t, - tf)e-^-'dt = e^-*°^/ (^) 

The spatial spectrum is given by the integral 

h=\ij{k)\^ - j \ij{k,uj)\^dujc^k-^ (26) 

For the surface elevations spectrum we obtain Phillips spectrum 




In our numerical experiments we observed the spectra both more steep for 
k —7- +00, and more shallow for k — oo than Eq. ( l26l) . So far, we have no 
proper explanation of this fact. 
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5. Turbulent quasibreathers in MMT model 

The Eqs. ©-([S]) have been solved numerically in periodic boundary con- 
ditions real space domain [0, 27r] for deep gravity surface waves case a = |, 
/3 = 3 and A = 1. Numerical integration has been performed through itera- 
tions of the implicit second order scheme in time and calculation of nonlinear 
term by Fast Fourier Transform technique. This numerical scheme preserves 
constants of motion of the approximated equation. 

To avoid high-frequency instabilities, the low-pass filtering has been ap- 
plied on every time-step through multiplication of the Fourier transform of 
the wave field by hyper-gaussian function, leaving about 90% of Fourier 
modes intact, while effectively suppressing the rest of potentially unstable 
high-frequency modes. Results were verified against the wave modes number 
change from 8192 to 16384 and 32768 for the same Cauchy problem. The cal- 
culations were continued typically up to thousands of the initial wave periods 
without loss of the accuracy. 

The initial condition was taken in the form of NLSE soliton 

^i^,0)= L \ (27) 
2k:^ cosh qx 



for km = 50, see Fig. HI It is known [l5| that simulation results essentially 



depend on the value of the nonlinearity parameter q/km- For q/k^, ^ 0.1 
the initial condition moves with the constant speed V without any notice- 
able shape change over characteristic length of at least dozens of simulation 
domain size 27r. For q/km > 0.1, the initial shape Eq.( l27|) starts to change in 
time and for q/km = 0.3 forms moving wedge-like growing structure with nar- 



rowing width. This behavior was interpreted in 15|] as possible collapse of the 
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initial condition over finite time, but further numerical simulation was not 
continued because of high-wavenumbers instability development in Fourier 
space, causing blow-up of the numerical scheme. In current research uti- 
lizing more sophisticated numerical approach, it was possible to follow the 
evolution of the same collapsing initial condition for practically unlimited 
time. We observed that, in fact, this collapsing initial condition evolves into 
localized non- stationary solution, periodically recurring to its initial shape. 
By analogy with cubical NLSE, it was interpreted as a breather-like struc- 
ture. 

The observed phenomenon is quite interesting: at q/km ~ 0.3 the ini- 
tial condition Eq. f l27|) evolves into localized object, but with "inner life". 
The shape of this object and the form of its spectra demonstrate irregular, 
stochastic behavior, which can be interpreted as some "intrinsic turbulence" . 
Time evolution of real space maximum of the solution is presented on Fig. 
O One should note that oscillations are quasi-periodic and their amplitude 
slowly diminishes in time, at least partially due to destruction of the breather 
by surrounding noise - that's the reason why we called this localized state 
by quasibreather. Almost identical picture of oscillations is seen from the 
second curve on FigJSl which presents the behavior of the second moment 
as a function of time. Both curves on Figl5] clearly indicate the presence of 
nonlinear oscillating structure in the wave system. 

FigJS] shows dependence of the frequency of these oscillations on the their 
mean level. The frequency has a tendency to grow with the growth of the 
oscillations level. This fact is in a correspondence with frequency dependence 
on nonlinear frequency shift. 
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FigJT] presents real and Fourier space of the system at t = 38.88, corre- 
sponding to the first maximum from FigJS] The real space picture of | (x) p 
shows that initial condition moved to the right with respect to initial condi- 
tion, growing in the amplitude and narrowing in width. Also, small portion 
of the initial condition has been separated in the form of the hump of much 
smaller amplitude. Fourier space contains two maxima: the right major peak 
approximately at km = 50, corresponding to the quasibreather, and the left 
smaller peak corresponding to the solution of the Eg. (120 p : 

ko = -(3 -V8)-km^ -8.6 

As shows Figd the spectrum remains localized near initial wave number 
k ~ km- This fact can be explained by conservation of both wave action and 
momentum. Thus, the turbulence inside the solution can be interpreted as 
an "envelope turbulence" . It is interesting that the area of this turbulence is 
localized both in real and Fourier spaces. 

Comparison with initial data shows that the spectrum gains power-like 
tails Ik — k~^'^ for negative k and Ik — k~^'^ for positive k. Recall that the 
simple collapse theory predicts Ik — k~^. Anyway, appearance of power- like 
tails indicates violation of smoothness of ilj{x,t). 

The observed singularity is of weak-collapse type. It is confirmed by the 
fact that the amount of Hamiltonian absorbed during 11 periods of oscil- 
lations of quasibreather (see Fig. E]) is approximately equal to 0.03% of its 
initial value. One should note that observed picture is universal: another 
snapshots of the system, taken at the times corresponding to subsequent 
maxima from Fig. [5l reveal the pictures similar to observed on Fig. [7] (see, 
for example. Fig. 
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FigJH] presents real and Fourier space of the system at t = 259.91, corre- 
sponding to the third trough from FiglSl The real-space picture of 
shows that amplitude of quasibreather has been diminished with respect to 
the state corresponding to Figd Fourier space exhibits both similarities and 
differences being compared to bottom of the FigjTl there are the same right 
main peak approximately at km = 50 and the left smaller peak approximately 
at ko = —8.6, but high-wavenumber tails decay much faster than power law. 
It means that ip{x,t) is smooth at the moment of minimum. 

For the illustration of the quasibreather temporal behavior, we present 
Fig. [TOl showing two states of the system taken at the moments when qua- 
sibreather reaches it's maximum and minimum amplitude in semi-log scale. 
It's quite obvious that spectral tails decay exponentially at the moment corre- 
sponding to the amplitude minimum of quasibreather, and decay as a power 
of wave number at the moment of the quasibreather amplitude maximum. 
This solution, therefore, periodically "breathes" between states of singularity 
formation and its regularization. 

Fig. [m presents surface elevation Eq. ([7]) for the same time as Fig. |H 
This picture looks qualitatively similar to experimentally observed "Three 
Sisters" killer wave on the ocean surface Ij and the resent results on freakon 
simulation on the deep water surface [9| . Fig. [12] shows slope elevations, 
corresponding to Fig. [11] These slope elevations values have the meaning of 
the original Euler equations for deep water surface gravity waves. 

One remarkable feature of the observed quasibreather is its co-existence 
with surrounding noise environment, associated with the radiation at the 
secondary spectral peak at fco = — 8.6. In fact, the surrounding weakly- 
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nonlinear noise could consists not only of radiation at wavenumber = —8.6, 
but also of the products of the initial condition decay into quasibreather and 
other waves. However, the wave action density in this noise is so small with 
respect to energy density in quasibreather, that this noise certainly cannot 
be interpreted as a kind of " condensate" . 

To analyze this situation, we performed the following experiment. In 
the middle of the simulation the real-space, containing quasibreather and 
surrounding noise, was "cleaned-up" through zeroing the function iIj{x) ev- 
erywhere except the carrier domain of the quasibreather. As a result, fur- 
ther evolution of the system starting from such "cleaned" initial conditions 
didn't show any qualitative difference from previous behavior - we observed 
immediate appearance of the surrounding noise at = —8.6 of the same 
characteristic amplitude, as we have seen before the "cleaning" of the real 
space. 

This observation lead us to the conjecture that quasisolitons and quasi- 
breathers exist only in quasi-equilibrium with weakly nonlinear wave noise 
environment. 

Another important observation, which distinguishes quasibreathers from 
oscillations of perturbed NLSE solitons, is periodical singularity formation 
at every time quasibreather reaches its maximum. This property is illustrated 
by both Fig. [7] (corresponds to the first maximum from Fig. [5]) and Fig. [9] 
(corresponds to the maximum number six from Fig. E]). 

In a nutshell, the gravity surface waves MMT model shows periodic fo- 
cusing of the initial condition Eq.f l27|) with weak-collapse singularity forma- 
tion exhibiting itself in power spectral tails and weakly nonlinear radiation 
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at secondary spectral maximum at ko = —8.6, which differs observed quasi- 
breather from previously known breather-like structures. The similarity of 
observed quasibreather in terms of water surface elevation with experimental 
"Three Sisters" wave packet and numerically observed freakon shows that 
even simplified model of gravity surface waves as MMT catches significant 
properties of the original exact equations. 



6. Conclusion 

On the base of numerical experiments we see that quasisolitons in frame 
of defocusing MMT model with parameters a = 1/2, /3 = 3 and A = 1 
are robust long-living objects, existing for hundreds of leading wave periods. 
Quasisolitons of large amplitude turn to quasibreathers. Their amplitude 
and spectral shape oscillate in time. These oscillations are accompanied by 
formation of weak collapses which can be compared with "white capping" of 
real ocean waves. 

We conclude that the "solitonic" scenario of freak waves is based on the 
equal foot with alternative " instantonic" scenario. We need to perform more 
numerical experiments in the frame of exact Euler equation to establish what 
scenario is closer to reality. 

Let us mention that oscillatory effects in solitons propagating on zero 



background were observed in paper 20|. However, in this paper the authors 



studied not single NLSE, but the system of coupled NLSE. The dynamics 
of this system is much more complicated. 
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8. Appendix I 

Now we address the following question - what value of A has to be chosen 
to provide the best possible modeling of real surface gravity waves on deep 
water? 

To answer this question, we notice that weakly nonlinear gravity waves on 
deep water surface with gravity acceleration g = I are described by so-called 
"Zakharov equation", which is exactly Eq.(jl]) at a = 1/2. 

The "real" coupling coefficient Tkkik2H ^ complicated homogeneous 
function of the third order: 

rpR 3rpR 

-'- ekeki€k2ek3 — ^ ^ kkik2k3 

ilxplicit expression for "real" Tj^k^^^j^^ was found, for instance, in the paper 
19|. 



n 

Functions T^^^j^g from [19] and Tkkikihi, given by Eq.(l5]), are essentially 
different. However, we can make them coincide in one point k = ki = k2 = 
by the proper choice of A. 
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According to 19 1 



R '^3 

Tkkikik-s = (29) 

But in the cited paper we used the "symmetric form" of the Fourier trans- 
form. If we define the Fourier transform according to Eq.Q, we must replace 
Eg. (1291) to 

'^kkkk = (30) 

Hence, to reach the best approximation to reahty, we have to put 

Tkkkk = (31) 

It means that we must choose A = 1. Then the shape of the surface ri{x,t) 
defined by Eq.Q is a model (rather approximate, of course) of a real water 
surface. From Fig. [11] one can conclude that the steepness of our breather is 
fairly high and hardly can be described by NLSE. 

9. Appendix II 

The vast majority of surface waves physical characteristics measurements 
is coming from stationary installations like oil platforms, presenting the water 
surface elevations as a time series. 

The water surface elevation itself is not a measure of the system non- 
linearity degree, since underlying equations are invariant with respect to 
stretching transformations, therefore surface waves of height varying by the 
order of magnitude can be of the same degree of nonlinearity. 
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The real physical characteristic of nonlinearity is the wave slope /i, which 
needs to be recovered from the surface elevations time series. Here we preset 
such simple estimate. 

By definition, the slope (same as steepness) is /i = fca, where k and 
a are the characteristic wave number and amplitude correspondingly. The 
connection between wave period T and wave number is 

, 47r2 



(32) 



For the 



amous "Draupner Wave" (also known as "New Year Wave", see 
FiglSl 2l|), T = 12 sec, a = 13.7 m and g = 9.81 m/sec, which gives fi ^ 0.38 
in accordance with our experiments. 
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Figure 1: Giant wave detected during a global census using three weeks of raw ERS- 
2 SAR imagette data, carried out by the German Aerospace Centre (DLR). This SAR 
data set was inverted to individual wave heights and investigated for individual wave 
height and steepness. The wave shown here has a height of 29.8 m. Adopted from 



http : / /www.esa.int/esaCP/ SEMOKQL26WD .indexA.html^^subhead'i 
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Draupner wave record 

January 1 1995 at 15:20, hs = 11.9m 
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Figure 2: Freak wave event detected from the Draupner oil platform on Jan.l, 1995. 



Adopted from littp://www. matli.uio.no/~karstent/seminarV05/Haver2004.pdf 
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Figure 3: Example of the situation when defocusing quasisolitons are possible. The disper- 
sion relation is a; = |/c|" for a < 1, f2 is negative and V is positive. The straight line always 
crosses the dispersion relation uj = uj{k) and, therefore, the denominator fl — kV + uj(k) 
in Eq. (jH) has zero. Quasisoliton takes place only in the defocusing case A = +1. 
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Figure 4: Real and Fourier space distributions of wave field. Top graph: \'tp{x,t)\'^ as a 
function of x fov t = 0. Bottom graph: Fourier spectrum logio as a function of 

signed logarithm of waves number sign(fc) log^Q |A;| for time t = 0. 
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Figure 5: Dependence of the solution maximum max {\ip {x, t))"^ , taken over integration 
domain [0,27r] (solid line, left axis) and the second moment f{k — ko)'^\ipk\'^dk (dotted 
line, right axis), on time t. The average wave-number is defined as fco = ^/\^^\2dk 
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Figure 6: Dependence of quasibreather maximum oscillations frequency on the mean level 
of these oscillations < \tl}{x,t)\^ >. 
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Figure 7: Same as FigHJ but for t ~ 38.88, corresponding to the 1st maximum from FiglH 
Left slope of the spectrum is approximated by function ~ k^'^-^ (dotted hne), right slope 
is approximated by function ~ k~^^^ (dashed line). 
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Figure 9: Same as Fig. |4l but for time t = 479.0 corresponding to 6th maximum from 
FiglSl Left slope of the spectrum is approximated by function ~ k^^-^ (dotted Hne), right 
slope is approximated by function ~ k~^^^ (dashed line). 
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Figure 10: Comparison of two spectra logj^Q |'(/'fe(t)P for time t = 259.91 (solid line, corre- 
sponds to the third trough on the FiglS]) and time t = 479.00 (dashed line, corresponds 
to the six' peak on the FigjS]), plotted as a function of wave-number k. This picture 
demonstrates that the spectral tails "breath" between exponential and power-like states. 
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Figure 11: Surface elevation 'q{x,t) as a function of real space coordinate x for time 
t = 479.00, corresponding to Fig. |9] 
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